Remarks about the thermostatistical description of the HMF model 
Part I: Equilibrium Thermodynamics 



L. Velazquez 1 and F. Guzman 2 

1 Departamento de Fisica, Universidad de Pinar del Rio, 
Marti 270, Esq. 27 de Noviembre, Pinar del Rio, Cuba. 
2 Departamento de Fisica Nuclear, Instituto Superior de Ciencias y Tecnologia Nucleares, 
Carlos III y Luaces, Plaza, La Habana, Cuba. 
(Dated: October 15, 2009) 

Abstract 

In this series of papers we shall carry out a reconsideration of the thermodynamical behavior of the called HMF model, a 
paradigmatic ferromagnetic toy model exhibiting many features of the real long-range interacting systems. This first work is 
devoted to perform the microcanonical description of this model system: the calculation of microcanonical entropy and some 
fundamental thermodynamic observables, the distribution and correlation functions, as well as the analysis of the thermody- 
namical stability and the relevant thermodynamic limit. 

PACS numbers: 05.20.Gg; 05.20.-y 



I. INTRODUCTION 

The present effort deals with the macroscopic descrip- 
tion of a toy model with a great conceptual signifi- 
cance: the called Hamiltonian Mean Field (HMF) model 
0, 0, [|. Despite the HMF model is enough amenable 
for allowing an accurate numerical and analytical char- 
acterization, it exhibits many features observed in more 
realistic long-range interacting systems such as: violent 
relaxation, persistence of metastable states, slow colli- 
sional relaxation, phase transition, anomalous diffusion, 

etc 0, i 1, 0, S 1 M ED, El El El El El- This 

reason explains why the present system is considered as 
a paradigmatic toy model in the understanding of the 
thermodynamical and dynamical properties of the real 
long-range interacting systems [l3| . 

The Hamiltonian of this model system is given by: 

N -. N N 

H " = E Yl^ + 2 9 E E t 1 - cos ^ - > (!) 

i—l i—l j — 1 

where I is the momentum of inertia and g the coupling 
constant. This is just an inertial version of the ferro- 
magnetic XY model [13], where the interaction is not 
restricted to first neighbors but is extended to all cou- 
ples of spins. Alternatively, this system could also be 
interpreted as a set of particles moving on a circle and 
interacting via a cosine binary potential. We shall refer 
mainly in the present work to the ferromagnetic interpre- 
tation. Besides the total energy E = -ff/v, a very impor- 
tant microscopic observable is the magnetization vector 

m = (j2k=i m t) l N wnere m i = (cosflj, sin6»j). 

General speaking, the macroscopic characterization of 
this model system by using the standard thermostatis- 
tical treatment does not significantly differs from the 
usual properties of the short-range ferromagnetic models. 
However, the understanding of its dynamical features is 



still an open pro blem attractin g m uch attention in the 
last years M El, 0, El, MM ■ 

The present work is just the first part of a series of 
papers devoted to perform a re-analysis of the thermody- 
namical behavior of the HMF model with the aim to ob- 
tain a better understanding of its remarkable dynamical 
behavior. As expected, we begin to address the equilib- 
rium thermodynamical properties of this model system 
by using the microcanonical description. As elsewhere 
discussed, the microcanonical ensemble is just a dynam- 
ical ensemble associated to the ergodic character of the 
many-body nonlinear Hamiltonian dynamics. Thus, the 
using of this statistical ensemble allows to keep a close 
relationship with the dynamical behavior of the HMF 
model studied by means of numerical microcanonical sim- 
ulations. Besides the entropy and some relevant thermo- 
dynamical parameters, the magnetic susceptibility, the 
particles distribution and the two-body correlation func- 
tion will be also obtained. Two novel aspects accounted 
for at the end of the present discussion are the analy- 
sis of thermodynamical stability and the relevance of the 
thermodynamic limit. 

II. EQUILIBRIUM THERMODYNAMICS 

The thermodynamical description of the HMF model 
by using the mean field approximation was considered 
in refs. [H, EH • We shall perform in this section the mi- 
crocanonical description directly working on the iV-body 
phase space. As usual, the microcanonical expectation 
value of any microscopic observable A (X) is obtained 
from the expression: 

A m = ^Sp[A(X)S(E-H N )] ) (2) 

where is microcanonical states density, f2 — 
Sp[S(E — Hn)], while Sp [Q] represents the N-body 
phase space integration: 



dXQ(X) = — J ^-^Q{9,L). 



(3) 



i N 9d N 
(2nH) 

where the factorial term considers the particles identity. 



A. Thermodynamical parameters 

The microcanonical states density is calculated as 
follows: 

1 f d N 9d N J 
n(E,N;I,g) = — J ——-S[E-H N (9,L-I,g)}. 



The integration by d L yields: 



(4) 



d 9 [E-V N (6;g)]^ 



\ N 



The accessible volume W is expressed in a dimensionless 
form as W = Slg/2. This last equation is rewritten when 
N tends to infinity as follows: 



W 



2ire 3 Ig 
NH 2 



lJV-1 



d 2 xn [m 2 + 2u - l] 2 /(m;7V), 

(6) 

where the dimensionless energy u = E/gN 2 and the mag- 
netization distribution function / (m;iV): 



f(m;N)=N 2 J 



d N 9 
(2,f 



fe=i 



(7) 



were introduced. The distribution function can be 
rephrased by using the Fourier representation of the 
Dirac delta function as follows: 



(2^) 2 



exp (iNK • m) I (-iVK 2 ^ 



N 



(8) 



where K = k + ix, being x a real bidimensional vector 
and /„ (z) the modified Bessel function of n-th order: 



In{z) 



2n 



+ 71 



exp (z cos 9) cos (nO) dd, (9) 



which satisfies the recurrence relations: 



2n. 



In+l (z) - I n -l (z) = In (z) , 



J„_i (z) + I n+ i (z) = 2 



dz 



(10) 
(11) 



The main contribution of the integral when N tends to 
infinity is obtained by using the steepest descend method: 



/ (m; N) ~ exp {— N [xm — In Iq (x 

2 



7^ ) ,; » '•''"''^ 



(12) 



where the functions k\ (x) and k,2 (x) are given by: 



ki (x) 



2m (x) 



+ — m (x) > 0, k 2 
x ax 



m (x) 



>0. (13) 



being m = |m| the magnetization modulus related to 
x = |x| by: 



m (x) 



h(x) 



(14) 



The validity of the above approximation is ensured by 
the positivity of the argument of the logarithmic term 
in Eq . (fT2|) . The derivation of this expression appears in 
the appendix IA II Notice that the integration function 
of Eq.© only depends on the modulus of m, and there- 
fore, the microcanonical expectation value of m identi- 
cally vanishes as a consequence of the nonexistence of 
a preferential direction of this vector. Nevertheless, the 
expectation value of m can differ from zero as a con- 
sequence of the occurrence of a spontaneous symmetry 
breaking (see in the next subsection). 

The main contribution to the entropy per particle 
s (u, N; I, g) = In W/N when N tends to infinity can be 
obtained by using again the steepest descend method as 
follows: 



s (u, N; I, g) — sq + max | — In [m 2 (x) + k] 

-xm{x)+]nI (x)} + O(— ), (15) 

where k = 2u — 1 and the additive constant so = 
i In (2we 3 Ig/Ny The stationary condition is given by: 



m (x) 



m 2 (x) — k 



m 2 (x) + k dx 



m(x)=0, (16) 



arriving in this way to the magnetization modulus versus 
energy dependence: 



m = mix) 1 • , m \ 

u = \m (x) fx + \ [1 - m 2 [x)\ ) ^ X G [ °' °° ) ' 
and m = if u > u c = 0.75. 

(17) 

The caloric curve is obtained from the canonical param- 
eter r) = gN(3 = ds (u, N; I, g) jdu: 
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FIG. 1: Microcanonical description of the HMF model in the 
thermodynamic limit: the caloric T (it) (solid line) and mag- 
netization m (u) (dash line) curves clearly reveal the existence 
of a second-order phase transition at u c — 0.75 with T c = 0.5. 
The divergence of the correlation function g (it) (dash-dot- 
dot line) manifests the existence of a long-range order at the 
critical point u c . 



_ 1 _ J x/m (x) 
T 



with < u < u c 



^ l/(2it— 1) otherwise, 



(18) 



where T is the dimensionless temperature. The depen- 
dences shown in FIG Q] confirm the existence of a continu- 
ous (second-order) phase transition at the critical energy 
u c = 0.75 with critical temperature T c = 0.5: from a fer- 
romagnetic phase when u < it c , towards a paramagnetic 
one with u > u r . 



Since m (x) drops to zero when u 
power expansion: 



the use of the 



t s 1 1 3 1 5 

to (x) — —x x H x + .. 

y ' 2 16 96 

allows to express x in terms of u as follows: 



32 . . 7 /32 

y K - M) + 3oly 



(19) 



(20) 



and therefore, the temperature and magnetization de- 
pendences in the neighborhood of the critical point with 
u < u c are given by: 



TM 1 2 ^64 a 



m (u) 



2 ( 4 

2\l -{u c ~ u) I 1 - — (ti c - «) + 



(dT/du) 1 undergoes a 



c(«+)-c(«-) 



(23) 



being C (u) = 1/2 when it > tt c . 

The results obtained so far are in a total agreement 
with the one derived from the mean field approxima- 
tion carried out in ref. 13]. As already evidenced, the 
thermodynamical features of the HMF model does not 
essentially differ from the ones exhibited by other ferro- 
magnetic systems in the mean field approximation [13]. 
Therefore, the solution above provided should lost its 
validity in the neighborhood of the critical point u c 
HH Ej EH- In order to show this fact, let us consider 
the second derivative of the maximization problem (| 1 5() 
which allows to check the stability of the stationary so- 
lution (fT7|l : 



h (x; it) = to' (x) 



to 2 (x) 



(to 2 (x) + k) 



1 to' (x) . (24) 



We have used the recurrence relations (TIT))) and (TTTj) 
in order to express the first derivative of the magne- 
tization as to' (x) = 1 — to (x) jx — m 2 (a;). The sta- 
tionary condition (fTB)) allows to substitute the relation 
k = k (x) = m (x) jx — m 2 (x) into (|24)l . Thus, the sec- 
ond derivative (|24[) is negative everywhere with the ex- 
ception of the stationary solution x = corresponding to 
the paramagnetic phase: 



h (x = 0; it) 



2u-l 



(25) 



which turns unstable when it < tt c . The function h (x; it) 
is related to the average square dispersion of the magne- 
tization throughout the expression: 



(<rm 2 ) ~ (to' (x) Sx) 2 = - 



(m'(x)f 
Nh (x;u)' 



(26) 



which diverges at the critical point of the second-order 
phase transition it c , yielding: 



(,5m 2 ) 



12it-l 



(27) 



when u > u c . Since the average magnetization square 
dispersion is directly related to the correlation functions 
gij — (m.j • m.j) — (mj) • (m^) via the formula: 



N(Sm 2 ) 



N ^ 



9ij = 9 (u) 



(28) 



the divergence of JSm. 2 ) implies the existence of a long- 

(21) range order at u c [13] (see also in FIGQ]). It is well-known 
that the existence of such a phenomenon significantly 

(22) modifies the behavior of the thermodynamical quantities 
close to the critical point u c . 
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B. Magnetic susceptibility 

The Hamiltonian of the HMF model is invariant under 
the translation operation 6^ — > 6k+ip, which is equivalent 
to the f7 (1) rotational symmetry acting on the rotator 
directions, whose existence implies the vanishing of the 
expectation value of the magnetization m. The U (1) 
symmetry is broken by modifying the Hamiltonian ([1]) 
with the incidence of an external magnetic field H as 
follows: 



H N -> H* N = H N - 9 NJ2 H ■ m*. 



(29) 



The vector H introduces now a preferential direction for 
the average magnetization, which leads to a modification 
of Eq. (J6]) as follows: 



W oc / d 2 m f (m; N) (to 2 + re + 2H • m 



AiV-i 



(30) 



as well as the following dependences of the caloric and 
magnetization curves: 
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FIG. 2: Microcanonical description of the HMF model un- 
der the incidence of a magnetic field: caloric (solid lines) and 
magnetization curves (dash-dot), as well as the microcanon- 
ical susceptibility (dot lines) for (a) with H = 0, (b) with 
H = 0.02, and (c) with H = 0.08. 



e(x;H) = 



TO (x) + H 1 r 9/nT , ,n ,n 

+ - [l - m 2 (x) - hm (x) , (31) 



2.T 



V=7P = 



T m (x) + H 



to (x) 



h{x) 
h (a) ' 



(32) 



where —H < e < +oo, < T < +oo, < m < 1 when 
x goes from the infinity to zero, being H = |H|. No- 
tice that we have distinguished between the dimension- 
less total energy e — u — Hm associated to the modified 
Hamiltonian (|29p and the dimensionless internal energy 
u of the Hamiltonian ([1]). These curves are represented 
in FIG [2] by considering different values of the external 
field H. Thus, the system exhibits a nonvanishing mag- 
netization m for every finite energy u, as well as there is 
now a smooth dependence between T and e. 

The caloric and magnetization curves (jTTJ) and (ITS)) 
are obtained from (|3ip and (|32p when H — > 0. The ex- 
pectation value of the magnetization drops to zero with 
the vanishing of H when u > u c , and therefore, the sys- 
tem recovers the original U (1) symmetry. However, a 
nonzero magnetization along the direction of the exter- 
nal field survives when the vanishing of H is carried out 
by keeping fixed its orientation. It means that the initial 
symmetry U (1) is spontaneously broken when u < u c . 

Let us now obtain the magnetization dependence for 
low values of H when e > u c . Since m (x) and x simulta- 
neously vanish when H — * 0, we are able to use the power 
expansion (|19p in order to obtain the power expansion of 
the H in terms of x starting from Eg. ([31]) : 



(33) 



H = 2x (e - u c ) + x 3 ^2e - j£\ + O (x 4 ) 
whose inversion leads to the expression: 



' H-l**^°H> + 0(H*). (34) 



2 (e - u c ) 



4 4 tr- 



ie - u c ) 



We finally obtain the m versus H dependence by substi- 
tuting this latter expansion into Eq.(|19p. yielding: 



m (e; H) = Al } - - ^ £ 2 \ 4 H 3 + O (H 4 ) 



4(e-e c ) 4 4 (e-u c ) 4 

(35) 

Thus, the microcanonical susceptibility x^™' when u > 
u c is given by: 

Am) _(dm\ _ 1 3 36e - 22 2 , 3 , 

(36) 

which diverges at the critical energy u c . 

The dependence of the magnetization m on the ex- 
ternal field H when e < u c can be represented as 
to = m(xn), where xh = x (e; H) is the solution of 
the problem: 



e [xh] H] = e < u c . 



(37) 



being e (x, H) the function defined by Eq. (|3"Tj) . The solu- 
tion for small values of H can be obtained by using the 
power expansions: 

x H = x + Cl (x)H + c 2 (x)H 2 + .., (38) 

m H = m (xh) = m (x) + xi (?) H + \2 (x) H 2 + .., 

(39) 
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being s(x, 0) = e and \i ( x ) = to'(x)ci(x), X2 (&) = 
|m" (x) (x)+to/ (x) C2 (x), and so on. The substitution 
of Eqs.([38f and ([39]) into (JSJ) yields the following results 
for the function XV 



Xi (x) 



m' (x) (m (x) - i) 



(-to (x) to' (x) - + ^m' (x)) 

and X2 = Pi (a?) M (%), being: 



/ \ / n . 1 2 / n, m{x) fxi{x) 

Pi (x) = Xi (x) + ~Xi O) 3- — TTT 

2 x a \m' (x) 



1 xi Qc) 

2x 2 to' (x) 



(Xi (x) + 1) 



1 1 m (x) 
p 2 (x) = -m(x) + — - — g . 

2a; 2x z to (x) 



(40) 



(41) 



(42) 



The function xi (x) represents the zero-order approx- 
imation of the power expansion for the microcanonical 
susceptibility Xh"^ i x , H) in terms of the magnetic field 

H. It is possible to show that Xh — Xi a l so diverges at 
u c as follows: 



It worth to remark that the microcanonical suscep- 
tibility Xh a l so admits negative values when e < ui- 
According to the well-known theorem derived from the 
canonical description: 



(c) 

Xh 



d(M) c 
dH 



(3G C 



(46) 



being G c = (AM 2 ) the average square dispersion of the 
magnetization within this ensemble, the magnetic suscep- 
tibility should be nonnegative. Actually, the microcanon- 
ical susceptibility Xh obtained in this subsection pro- 
vides a measure of the magnetic sensibility of the system 
at constant energy instead of at constant temperature. 
It can be shown that the microcanonical counterpart of 
the identity (I46[) in the thermodynamic limit plj is given 
by: 



(m) 

Xh 



d(M) T 
dH 



= f3G m + (M) r 



d(M), 
dE 



(47) 



which clarifies that the negative values of Xh^ come from 
the term MdM/dE since usually DM /dE < in the 
ferromagnetic phase. 



Distribution functions 



(m) 

Xh 



1 1 



8 (ii c - e) ' 



(43) 



when e < u c . 

The above estimations of the microcanonical suscep- 
tibility Xh^ are onr y applicable when 1x1-^1 << !■ I n 
general, the magnetic susceptibility can be obtained from 
the formula: 



(m) 

Xh 



to' (x) 



( dx 
(dH 



(44) 



Let us now obtain the microcanonical n-body distribu- 
tion functions of this model: Fm^ [y] = (yi, ■ ■ ■ Vn)'- 

M = ^Sp [S ( yi - Xl )...5(y n - x n ) S (E - H N )} , 

(48) 

where Xk = (Ok, Lk). We shall assume that in the fer- 
romagnetic phase there is a nonvanishing magnetization 
m = (m, 0), which could be obtained by considering the 
Hamiltonian (j2"9"| with H= (H, 0) and sending H to zero. 

The state density fi can be rephrased in a functional 
form T as follows: 



which yields: 



(m) 

Xh 



1 

2x 



\2x 



H 



2x 2 



(to + H) to' 



(45) 

The microcanonical susceptibility Xb^ was a l so re P re " 
sented in FIG[2]for different values of the external field 
H . Notice that Xh rema i ns almost constant at negative 
value Xh 1 '' — — 0.46 when e < u\ — 0.5. The suscepti- 
bility begins to grow beyond the point u\ and vanishes 
at ui ~ 0.61 independently of the value of the external 
field H , a behavior provoked by the presence of the factor 
m(x) — 0.5/x. The system exhibits a large susceptibil- 
ity under the influence of a very small external magnetic 
field in the energetic range (u2,u c ), which is related to 
the existence of the second-order phase transition at u c . 



f2 oc.F[^,iV] = J <fm f (m,N) [tp]^ 



(49) 



where ip — tp (m; N) = 2u — 1 + m 2 (see in Eq.©). It 
is easy to see that the n-body distribution functions can 
be also rephrased by using this same functional form as 
follows: 



where 5tj) n is given by: 



{n) ^ T[^ + 5ip n -N -n] 



(50) 



N 
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being m fc = m (9 k 
vector, 5mk = mfc 
and it — m 
k-th rotator. 



(cos 8k, sin Ok) the magnetization 
m, the corresponding dispersion, 
= \p\ — m ■ rrife = £fe the energy of the 
We also consider for convenience the di- 



mensionless momentum p k = Lk/Lo by introducing the 
characteristic unit Lq = y/IgN . 

The expression (j50|) suggests to perform a perturba- 
tive expansion in power series of 1/N of the functional 
In T (ip + 5i/) n ; N — n) by taking into account the Gaus- 
sian localization of these integrals when N tends to in- 
finity. We obtain after some algebra the following result: 



d 2 m f(m,N- n) ^("-"J-^W [ y - m ] , 

(52) 




-V 2 (E e fe) +(n + 2)vJ2 e >< 

k=l 



\k=l 



0I N 



(53) 



being r\ — ip , the dimensionless inverse temperature. 
The progressive calculation demands a refinement of the 
steepest descend method in order to account for the 1 /TV- 
contributions obtained beyond of the Gaussian estima- 
tion (see in appendix I A 2|) . 

The zero-order approximation of the one-body distri- 
bution function is given by: 



f m (9, p; u) = C cxp (-tie) , (54) 

where e = e(6,p) = \p 2 — m • m (9), the normalization 
constant is given by: 



C~ L = \ / — 2ttIq Onm) . 

Tj 



(55) 



and the magnetization vector m satisfied the self- 
consistent relation: 



m = J d9dp m (9) f m (9, p; u) 



(56) 



FIGE] shows the angular distribution function p (9; u) for 
different energies: 



p(9;u)= dp f m (9,p;u) = , (57) 



2tt/ (x) 



where the relation x = ?ym was taken into account. No- 
tice that the ferromagnetic states with u < u c are char- 
acterized by the existence of a clustered distribution of 




FIG. 3: Angular distribution function p (8; u) for different 
energies values. The clustered distributions in the ferromag- 
netic phase progressively become in a uniform distribution 
when the energy u > u c . 



the angular variables 9k around the direction of the mag- 
netization vector m, while a uniform distribution is ob- 
served in the paramagnetic phase with u > u c . This is 
the reason why the ferromagnetic states are referred as a 
clustered phase, while the paramagnetic ones are referred 
as a homogeneous phase. 

The terms of the two-body correlation function 

gin (xi,X2) = Fm' (xi,x 2 ) - (%i) Fm (x 2 ) surviv- 
ing the first-order approximation are given by: 



1$ (xi,x 2 ) 



Jjhh [$12 



rfg (u) 6 



12 



(58) 

where fa = f m (6i,Pi;u) is the zero-order approxima- 
tion of the one-body distribution function (l54l) . g (u) 
Y (Am 2 ), the correlation function derived from the av- 
erage magnetization dispersion (Am 2 ) showed in FIGHJ 
and $12 and 612, two new form factors given by: 



$12 = rjSmi ■ Sm 2 — 2rj 2 5eide 2 , 
©12 = [2rjmSei + Smi) ■ (2rjmSe 2 



Sm 2 ) 



(59) 
(60) 



The quantities Sei = e% — (e) and Jnij = m^ — m rep- 
resent the energy and magnetization desviations respec- 
tively. The presence of the function g (u) in Eq. (f58|) leads 
to the divergence of the two-body correlation function 
at the critical point u c . The spacial two-body corre- 
lation function in the homogeneous phase is given by 
9 ($i> $2) = (2tt)~ 2 c 2 (0i, 9 2 ;u), where the function: 



)= 1 + 7 



1 1 



4 u 



1 



2u- 1 



■ cos ( 



h-0 2 ), (61) 
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diverges at the critical point in terms of the inverse tem- 
perature 77 as follows ~ 4 (rjc — 77) _1 cos (9% —62)- This 
asymptotic behavior is consistent with the one estimated 

in the refs. mm 13, but flUD is now the exact micro- 
canonical result within the first-order approximation. 

The inexistence of three-body terms in Eq. ((53|) 
straightforwardly leads to the vanishing of the three-body 
correlation function 5™' (xi, x 2 , x 3 ): 



( 3 ) I \ c( 3 ) I \ ( 2 ) / ^ Z?(!) I \ 

gin' {Xl,X 2 ,X 3 ) = Frn' (Xi,X2,X 3 ) -gin {Xl,X 2 )Fih (x 3 ) 

-g$ (x 2 ,x 3 )F^ (x^-gtt (x a)Xl )F<P (x 2 ), (62) 

in the first-order approximation, and hence, the three- 
body correlations are just 0(l/N) size effects. Such 
a microcanonical result will be taken into consideration 
during the derivation of suitable dynamical equations for 
the distribution and correlation functions based on the 
well-known BBGKY hierarchy. 



III. THERMODYNAMIC STABILITY 

Thermodynamic stability concerns to the question 
about when a given admissible macrostate characterized 
by a certain energy and magnetization is stable or un- 
stable under thermal fluctuations associated to the ther- 
modynamic equilibrium where the temperature and the 
magnetic field act as constant control parameters. Such 
a question is also intimately related to the nature of the 
correspondence between the control parameters (77, H ) 
and the controlled system observables (u, to) (internal en- 
ergy u and the projection of the magnetization to along 
the magnetic field), that is, the existence or inexistence 
of the ensemble equivalence. 

The analysis starts from the consideration of the parti- 
tion function Z = J dXpj exp [— (3H^\ and the introduc- 
tion of the Planck thermodynamical potential per parti- 
cle p(n,H) = —\nZ(i],H) /N, which can be rephrased 
as follows: 

exp [— Np (n, H)] oc J dmduexp{—Np(r],H;u,m)}, 

(63) 

where p (77, H; u, to) = r]u + Ato — s (it, to) with A = —[3H 
and the detailed entropy s (u, to) given in a parametric 
form by: 

s (u, m) — s + - In {k + to 2 (x) } — xm (x) + In J (x) 

(64) 

within the approximation provided by the steepest de- 
scend method. The exponential function of the integral 
(|63|) exhibits sharp peaks around its maxima when N is 
large enough, a behavior that allows to rephrase the inte- 
gral (|63[) within the Gaussian approximation as follows: 



y^exp {-Np k ) J exp 



NB (u, m; Uk, rrik) 



dmdu, 
(65) 



where p^ = p(n, H;uk,rrik) and the bilinear form 
B (u,m; u k ,m k ) = H uu Aul + (H urn + H mu ) Au k Am k + 
H m mAm\ with Au k = u — Uk and Arrik = m — rrik ob- 
tained from the entropy Hessian Hij : 



j r f Fl uu H mu 

lj ~ H H 

±J -um 11 rnrn 



d 2 s/du 2 d 2 s/dmdu 
d 2 s/dudm d 2 s/dm 2 



evaluated at the fc-th stable stationary point: 



ds (u k ,m k ) 
du 



A = 



ds (u k ,m k ) 
dm 



(66) 



(67) 



which satisfies the negative definition of the Hessian ma- 
trix (p)B"|) . From the thermodynamical point of view, such 
maxima represent the coexisting macrostates or phases 
appearing for given values of (3 and H . Obviously, there 
is ensemble inequivalence when there exist only one maxi- 
mum representing a unique stable phase with given values 
of (u, m). The greatest peak during the phase coexistence 
corresponds to the stable phase, while the other repre- 
sent metastable states. The greatest peak corresponding 
to the stable phase provides the main contribution of the 
integral (f63|) allowing to estimate the Planck thermody- 
namic potential per particle as follows: 



p (77, H) = inf {77M + Am — s (it, m)} . 



(68) 



This is just the Legendre transformation which consti- 
tutes a fundamental stone of the thermodynamic formal- 
ism [20(. All those admissible macrostates of the system 
satisfying the stationary condition (l67|) but do not obey 
the negative definition of the entropy Hessian matrix (|66|) 
are precisely the unstable macrostates. 

Performing the calculations, the canonical parameters: 



ds 



A 



ds 



du k + m 2 (x) ' dm k + m 2 (x) 



-x, (69) 



allow to express the magnetic field H and the tempera- 
ture T — /3 _1 in terms of the macroscopic observables k 
and m: 

H = H (k, to) = x (k + to 2 (x)) — to (x) , 

T = T(K,m) = K + m 2 (x) . (70) 



According to the entropy per particle (|64|) , the physically 
admissible regions satisfy the inequalities k + m > and 
to 2 < f. The regions of stability in the plane (k, m) are 
determined from the entropy Hessian: 

2 / -2 -2m 
\ —2m k — to — a [K + m ) 

being a -1 = to' (x) = 1 — to (x) jx — to 2 (x). The deter- 
minant D = dct Hij = — 2 (k + to 2 ) [l — a [K + m 2 )] 
vanishes at the boundary of the unstable region, a curve 
which is parametrically represented as: 



[k (x) , to (x)] — 



1 _ m(x) _ 2m2 (a . )jTO(x) 



(71) 
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FIG. 4: Magnetization curves of the HMF model and their 
stability in the plane (ft, m). Solid lines: magnetization curves 
at constant magnetic field H 7^ 0; Thick solid lines: magneti- 
zation curves at H = 0; Dashed lines: isothermal magnetiza- 
tion curves. Zones outside the curvilinear region PCQA are 
canonically stable. The dark zone inside this region is canon- 
ically unstable, that is, the region of ensemble inequivalence, 
while the white zones enclose the metastable states. Point C: 
critical point of the continuous phase transition. 
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FIG. 5: Caloric curves of the HMF model in the neighborhood 
of the critical point C. Solid lines: caloric curves at constant 
magnetic field H 7^ 0. Thick solid lines: caloric curves at 
H = 0, where BC-(PQ) and CA are the stable and unstable 
branches respectively. The nonanalyticity of the canonical 
thermodynamic potential at the critical point C follows from 
the bifurcation undergone by the caloric curve at H = in 
this point. The dash-dotted line divided the caloric curves 
into metastable and unstable regions. 



These results are summarized in FIG|4]and FIGiH The 
curve PAQ is the boundary of the physically admissible 
macrostates. The thick solid lines represent the magne- 
tization curves at H = 0. The stable branch PCQ con- 
stitutes the boundary between the stable and metastable 
regions. The critical point of the continuous phase tran- 
sitions C: (k c , 0) with k c — 2u c — 1 = 0.5 is just a bi- 
furcation point of the magnetization curves located at 
the endpoint of the metastable region which touches also 
the stable and the unstable regions, and therefore, it is a 
point of marginal stability. 

The unstable region (dark zones in the plane (k, m) 
shown in FIGf4|) is characterized by the presence of 
negative values of the magnetic susceptibility at con- 
stant temperature Xh = (dM/dH) T , a behavior canon- 
ically anomalous for this kind of model system in terms 
of the well-known thermodynamical identity Xh ~ 
/3(<5Af 2 ), but which is microcanonically admissible. Ac- 
tually, it possesses the same anomalous character of the 
macrostates with a negative heat capacity C = dE/dT < 
observed in several systems 22[ in terms of the ther- 
modynamical identity C = 1 (<5i? 2 ). 

The lost of analyticity in the thermodynamic limit of 
the thermodynamic potentials like the microcanonical en- 
tropy per particle s (e, H) = sup m {s [u(= e + Hm) , m}} 
or the Helmholtz Free energy per particle / (/3, H) = 
inf E {e — Ts (e, H)} = mf u ^ m {u — Hm — Ts(u,m)} at 
the critical point C can be related to the bifurcation 
of the magnetization curve shown in FIG0] or caloric 



curve show in FIG (5] for H = 0. The system is un- 
able to follows the "trajectory" BCA with m = since 
the branch CA is located inside the region of ensemble 
inequivalence where there exist anomalous macrostates 
with x'h 0. The large thermodynamic fluctuations 
existing there provoke a sudden change of the original 
tendency BCA following in this way anyone of the sta- 
ble symmetric branches CP or CQ with a nonvanishing 
magnetization m 7^ 0, and consequently, the occurrence 
of a spontaneous symmetry breaking [2fJ. 



IV. SOME FINAL REMARKS 

As already illustrated, the thermodynamic properties 
of the HMF model do not essentially differ from the other 
ferromagnetic models with short-range interactions. We 
have previously shown that the relevant microcanonical 
thermodynamic variable is u = E/gN 2 , and the charac- 
teristic energy is Eg = gN 2 . The thermodynamic limit 
is carried out when N is sent to the infinity by keep- 
ing fixed the dimensionless energy u. This is the same 
thermodynamic limit introduced in ref. [13j| in order to 
perform the mean field description of this model. How- 
ever, our analysis reveals that the entropy per particles 
is ill-defined in the thermodynamic limit: while the term 
of s (u, N; I, g) containing the relevant thermodynamical 
variable u is iV-independent, the A^-dependent additive 
constant sq = j In (2i:e 2 Ig/N^ diverges when N — > 00 
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(see in ref. [23] ) . It means that the ill-behavior of the en- 
tropy per particle is unavoidable without considering an 
appropriate iV-dependence for the coupling constant g . 

The using of an appropriate A-dependence in the cou- 
pling constant g is usually identified with the Kac pre- 
scription [24| . The standard usage of this procedure is 
to consider certain dependence g (N) that ensures the 
extensive growing of the total energy E, since the en- 
ergy per particle e = E/N is keep fixed in the thermo- 
dynamic limit N — > oo. This condition demands that 
Eq/N = gN = 7 = const. Although most of works 
devoted to the HMF model make use of this condition 
[I S S, 0, S, B M, ED, El d EI d , it is very easy to 
verify that its application does not avoid the divergence 
of the entropy per particle in the thermodynamic limit, 
limjv^oo so = oo. 

A closer look to this question clarifies that this proce- 
dure should be applied with care. Firstly, the coupling 
constants determine the characteristic temporal scales 
acting on the dynamical evolution of a given system, and 
therefore, any scaling A-dependence of these coupling 
constant affects the system dynamical behavior during 
the imposition of the thermodynamic limit N — > oo. This 
obvious remark is very important to take into considera- 
tion since the results of many numerical simulations evi- 
dence the noncommutativity of the thermodynamic limit 
N — > oo with the infinite time limit T — > oo necessary 
for the equilibration of temporal averages @, 0, H, Q : 

lim km (A) (T, N) + lim km (A) (T, N) , (72) 

(where (A)(T,N) = £ A[X N (t)]dt/T) as a conse- 
quence of the divergence of the relaxation time r eq in 
the thermodynamic limit, lirn/v^oo r eq = oo. In the au- 
thors opinion, the origin of this anomaly could be related 
to an inappropriate use of the Kac prescription. This ar- 
gument follows from questioning the main motivation of 
introducing the Kac prescription: to deal with an exten- 
sive energy. The extensivity of the energy is a thermody- 
namic feature of the extensive systems intimately related 
to the statistical independence or separability of a large 
system in independent subsystems appearing as a conse- 
quence of the incidence of short-range forces. Obviously, 
such a microscopic picture is outside the context within 
the HMF model, where the long-range character of the 
microscopic interactions implies its intrinsic nonexten- 
sive nature. 

The HMF model could be consider as a limit case with 
a — > of a parametric family of ferromagnetic models on 
a square lattice [25| whose potential energy is given by: 



(cos (6i — 6j)) w ) by using a lattice with spacing a and 
linear dimension R oc y/N : 

, : R 2-Krdr 

D a = V a /N ~ a 



= 27r /(a-2)-V--#-«) a ^2 

\ ln(JZ/a) a- 2 { ' ' 

allows to understand that such as models are extensive in 
the thermodynamic limit when a > 2 and nonextensive 
elsewhere. As a consequence of the existence of long- 
range correlations, the nonextensive system cannot be 
trivially divided in independent subsystems, and there- 
fore, there is no physical reason to justify the imposition 
of the extensive character of the total energy in this case. 

The using a suitable scaling dependence for the cou- 
pling constants in order to regularize the thermody- 
namic parameters and potentials of a given system in 
the thermodynamic limit might be applicable in some 
way that also ensures a well-defined dynamical behavior 
characterized by the commutativity of the limits limjv^oo 
and liniT^oo. Generally speaking, all that it is neces- 
sary to demand is the non divergence of the relaxation 
timescale t^, when N — > oo. With some recent excep- 
tions [ill, EH, most of numerical studies of the micro- 
scopic dynamics of the HMF model revealed a charac- 
teristic relaxation timescale for the HMF model given by 
T eq = T m i c N, being T m j c = y/I/Ng the characteristic mi- 
croscopic time of the rotators evolution [E @, S H, B EH • 
It is easy to realize that the imposition of the scaling 
g/N = 7 = const ensures both the well-defined behav- 
ior of the additive constant of the entropy per parti- 
cle So = \ m (27re 2 /<7/A r ) and the relaxation timescale 
r eq = yJIN I g when N — > oo. This ansatz leads to a 
power-low growing of the characteristic energy Eq oc A^ 3 , 
which does not involves in principle any physical incon- 
sistency since the HMF model is nonextensive. We shall 
return to this discussion in the forthcoming papers. 



APPENDIX A: DEMONSTRATIONS 
1. Derivation of Eq. (Tl2j) 

We perform the power expansion up to the second or- 
der approximation of k for the complex number z: 



2x 



i=l j= 



,y><-> 

--1 i] 



(73) 



1 (n-k) 2 + 0(|k| 2 ), 
(Al) 

being x = |x| and n = x/x. Let us now carry out the 
power expansion for the function ln/o (x + Az) in term 
of the variable Az up to the second order approximation, 
being Az = z — x: 



where r« is the lattice distance between the i-th and j-th 
rotator. The integral estimation of potential energy per 
particle in the paramagnetic phase (where the averages 



, t i . . , . . 1 dm (x) , . , 2 ^ 

ln/ (x + Az) = +m (x) Az+- ^ ' (Azy+O 



(Az; 

(A2) 
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where m (x) = I± (x) /Iq (x). Thus, the exponential func- 
tion: 



Eq. (|A7p allows to express the exponential function as 
follows: 



exp[7VK-m + lnio(z)], 



(A3) 



can be rewritten by dismissing the terms O ( |k| ) as fol- 
lows: 



exp (— N [xm — ln/ (x)]} x 



x exp 



N 



K\ (x) k\ — K2 (x) k\ 



(A4) 



where x || m and m = |m|; k\ = (n • k) and k\ = 
k 2 — (n ■ k) ; while k\ (x) = 2m (x) jx + dm (x) /dx and 
K2 (x) = m (x) jx. Since m (x) is an odd monotonic in- 
creasing function, with m (x) > when x > 0, the func- 
tions Ki (x) and K2 (x) are always nonnegative. The inte- 
gration of the expression (|A4[) leads to the result shown 
in Eq.Qr 



2. Refinement of the steepest descend method 

The calculation of the first-order approximation of the 
correlation functions involves a little refinement of the 
steepest descend method. The aim is the determina- 
tion of the l/A^-contributions of the average of a physical 
quantity A (x) as follows: 

(A) O = J A (x) exp [Ns (x; a)} (fx, (A5) 

where the partition function f2 is given by: 

Q = J exp [Ns (x; a)} cPx. (A6) 

Firstly, we perform the power expansion of the function 
s (x) around its maximum point Xo as follows: 

s (x) ~ s (xo) + c 2 ■ (Ax) 2 + R (Ax; x ) , (A7) 

where the higher-order contributions i?(Ax;xo) are de- 
noted by: 



i?(Ax;x ) = J2 c ™ ■ (Ax) r 



(A8) 



m— 3 



Hereafter, we shall consider the following convention in 
order to simplify the notation during the calculation. For 
example, the term c m ■ (Ax) represents the tensorial 
product: 



•(Ax) m = i-y — 



d m s (x ) 



m! ^ dx h dx i2 . . . dx l * 



-Ax 11 Ax 42 ...Ax 1 



(A9) 



exp [Ns (x)] = exp [Ns (x )] exp 



ic(Ax) 2 



F(Ax;x ) 



(A10) 

where uj = —2Nc2 and the function F(Ax;xq) = 
exp [NR (Ax; xo)] is given by: 

oo 

F(Ax;x ) = ^-iV fe i? fe (Ax;xo). (All) 



k=0 

This latter result can be conveniently rewritten as a 
power series expansion: 



f (Ax;x ) = K m - (Ax)' 



m=0 

where 

{]} \ k=3 ) k=3 J 

being d (n) an integer function defined by: 

1 n = Q 



(A12) 



(A13) 



6(n) = 



n^O 



(A14) 



It can be checked that the first 10 coefficients of the power 
expansion (|A12|) are given by: 

KQ = 1, Kl = K2 = 0, K3 = Nc 3l K4 = Nd, 

k 5 = Nc 5 , k 6 = Nc 6 + ^ 2 c 2 , k 7 = Nc 7 + N 2 c 3 c 4 , 
k 8 = Nc 8 + N 2 ^c 3 c 5 + ic 2 j , (A15) 
k 9 = Nee + N 2 (c 3 c 6 + c 4 c 5 ) + ^N 3 4, 
kio = Nc 10 + N 2 (c 3 c 7 + c 4 c 6 + ^yc 5 c 5 j + ^A f3 c 2 c 4 . 



Taking into account that in the Gaussian integration: 

1 



((Ax)") = N J d"xexp 



-uj ■ Ax 2 



(Ax) 



(A16) 

survive only the odd dispersions, denoting by <?2fc = 
N k ((Ax) 2 *), being Nq = wdet(-^w) the normaliza- 



tion constant, the first-order approximation of the parti- 
tion function f2 is given by: 



( - } - = \ 1 + (cm + ^c 2 3 g 6 



(A17) 



where Sl = exp [Ns (x )] /N is the zero-order approxi- 
mation of the partition function. 
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The calculation of the integral (|A5|) is carried out by 
performing the power expansion of the quantity A (x) 

around the maximum point x ^ q _ q q Z^ ) + — (a^ g 2 + 

oo 

A (x) = J2 A(m) ■ ( Ax ) m . (A18) + (A®* + A^ca) .94 + } , (A22) 



1 i9 m A ( x ) 1 



being: which leads to the following result 

d m A (x ) 

W 

. (A19) Restoring now the ordinary notation, we finally obtain: 
which allows to rewrite: 

A(x)F(Ax;x ) = E£n(Ax) n , (A20) (A) = A(x ) + ^gg(A^A^)+ (A24) 

being: 

A r > +- V ^ (xo) 93s(xo) AT (Ax* 1 Ax i2 Ax ia A3;' 4 ^ 

B»=V4 W ^ m . (A21) +3! 2^ ftc*! dx*dx*dx* { ' ' 

* ' il 12*3*4 

rn— 1 

Thus, the first-order approximation of the integral where (Ax^Aa^ 2 ) cx 1/jV and (Ax 11 Air 12 Ax l3 Ax u ) cx 
(|A5]) is given by: 1/N 2 . 
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